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Abstract 



The results of a search for W boson pair production in pp collisions at -y/i = 
1.8 TeV with subsequent decay to en, ee, and fifi channels arc presented. 
Five candidate events are observed with an expected background of 3.1 it 0.4 
events for an integrated luminosity of approximately 97 pb~^. Limits on the 
anomalous couplings are obtained from a maximum likelihood fit of the Et 
spectra of the leptons in the candidate events. Assuming identical WWj and 
WWZ couplings, the 95% C.L. limits are -0.62 < Ak < 0.77 (A = 0) and 
-0.53 < A < 0.56 (Ak = 0) for a form factor scale A = 1.5 TeV. 

PACS numbers: 14.70.Fm 13.40.Em 13.40.Gp 13.85.Rm 



3 



Typeset using REVT^]X 



Gauge boson self- interactions are a direct consequence of the non-Abelian SU{2) x U{1) 
gauge symmetry of the standard model (SM). The trilinear gauge boson coupling strengths 
can be measured directly by studying gauge boson pair production. Hadron collider ex- 
periments have established the electroweak coupling of the W boson to the photon |l| and 
the existence of the coupling between the W boson and the Z boson [0,0, and have placed 
constraints on anomalous WlVj and WW Z couplings . Measurements of the couplings 
also have been reported by the LEP collaborations H]. 



The WW'j and WW Z vertices can be described by a general effective Lagrangian |T0 



with two overall coupling constants g-ww-^ = ~e and gwwz = — e ■ cot 6w (where e is the 
W^ charge and 9w is the weak mixing angle) and six dimensionless coupling parameters 
gX , nv, and Ay (V^ = 7 or Z), after imposing C, P, and CP invariance. Electromagnetic 
gauge invariance requires that gj = 1. The effective Lagrangian becomes that of the SM 
when gl = gf = 1, i^v = = — 1 = 0), and Xy = 0. In order to preserve 

unitarity at high energies, the anomalous couplings are modified by form factors with a 
scale A (e.g. Xv{s) = Ay(0)/(1 + s/K^Y). Limits on these coupling strengths are usually 
obtained under the assumption that the WW'j and WWZ couplings are equal {gj = gf = 1, 
Ak^ = Ak^, and = A^), leaving two independent couplings to be determined. In 
another approach |TT| , the anomalous couplings are formulated in a framework that explicitly 
respects the SU{2) x U{1) gauge invariance, but contains more general terms than those of 
Ref. Comparison of the two formalisms leads to simple equations (HISZ relations) that 
relate anomalous couplings in the general effective Lagrangian. 

In this paper we present the results of a search for pp — >■ WW + X ^ ifuu' + X at 
■s/s = 1.8 TeV, where i,£' = e or fj,. Limits on anomalous WWj and WWZ couplings 
are obtained for both the equal couplings and the HISZ relations by a maximum likelihood 
fit of the observed two-dimensional spectra of lepton transverse energy E^- This method 
provides tighter limits on anomalous couplings than those from the measurement of the cross 
section The WW —>■ ii'uu' channel has significantly less background, albeit with a 

smaller branching ratio, and is more sensitive to WW production with the SM couplings 
than the WW/WZ — >■ ii'jj/iijj channel. Therefore, limits obtained from this analysis are 
complementary to those from the WW/WZ — > ii^jj/iijj analyses 

The data sample corresponds to an integrated luminosity of approximately 97 pb^^ 
collected with the D0 detector during the 1992-93 and 1993-1995 Tevatron collider runs at 
Fermilab. The results based on the 1992-1993 data sample of approximately 14 pb~^ were 
previously reported This paper describes the analysis of the 1993-1995 data sample 

and gives the combined results from the two analyses. 

The D0 detector consists of three major components: the calorimeter, tracking, and 
muon systems. A hermetic, compensating, uranium-liquid argon sampling calorimeter with 
fine transverse and longitudinal segmentation in projective towers measures energy out to 
|?7| ~ 4.0, where r] is the pseudorapidity. The central and forward drift chambers are used to 
identify charged tracks for \ri\ < 3.2. There is no central magnetic field. Muons in the central 
region are identified and their momenta measured with three layers of proportional drift 
tubes (PDT's), one inside and two outside of magnetized iron toroids, providing coverage 
for |?7| < 1.7. In addition, scintillation counters mounted on the outer layer of PDT's provide 
time information for muon identification and cosmic ray rejection. 

Event samples are obtained from triggers with the signature of leptonic W boson de- 
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cays. The ee and e/x samples are selected from events passing a trigger which requires an 
electromagnetic cluster with Ex > 20 GeV and missing transverse energy I^t > 15 GeV. 
The integrated luminosity for this sample is 82.3 ± 4.4 pb~^. The /i/x sample is selected 
from events passing a trigger which requires at least one muon track in the fiducial region 
of \ri\ < 1.0 with pt > 15 GeV/c, energy deposition in the calorimeter consistent with the 
passage of a muon, and no hits due to cosmic ray muons in the scintillator located outside 
of the muon chambers. The integrated luminosity for this sample is 65.2 ± 3.5 pb~^. 

Isolated electrons are identified using a likelihood function formed from four variables: 
the electromagnetic energy fraction of the calorimeter cluster, the of longitudinal and 
transverse shower shapes compared to test-beam and Monte Carlo electrons, the ionization 
energy deposition (dE/dx) in the central tracking detector associated with the matching 
track, and the distance between the projected track position and the centroid of the en- 
ergy cluster at the calorimeter. This likelihood function is used to discriminate between 
electrons and photon conversions, photon showers overlapped with a charged hadron track, 
and hadronic showers with large electromagnetic content. For a given identification effi- 
ciency, this method provides a background rejection 2-3 times higher than a method that 
places requirements on individual variables [|1^]. In the central region, |?7| < 1.1, the electron 
detection efficiency is (59.9±0.8)%; in the forward region, 1.5 < \rj\ < 2.5, it is (47.1±1.4)%. 

Muons are required to have associated hits in at least two of the three layers of the 
muon system. They must be isolated from jets (A7^(/i, jet) > 0.5 for E^r^^ > lOGeV, where 
A7?.(/i, jet) is the separation between muon and jet in rj — cp space) and have energy deposition 
in the calorimeter consistent with a minimum ionizing particle. The muon track is required 
to point to the primary event vertex within 25 cm in the plane transverse to the beam and to 
occur at a time, as measured by the PDT's, within 200 ns of the beam crossing. The muon 
detection efficiency is (70.1 ± 3.1)% within the fiducial acceptance of \ri\ < 1.0 employed in 
this analysis. 

The ee candidate events are selected by requiring the leading electron to have Et > 25 
GeV and a second electron to have Et > 20 GeV. The I^t is required to exceed 25 GeV. 
Background from events with high-p^ Z bosons is reduced by removing events with dielectron 
invariant mass within 15 GeV/c^ of the Z boson mass. The e/i candidate events are selected 
by requiring an electron with E^ > 25 GeV and a muon with px > 15 GeV/c. We require 
> 25 GeV (where ^f?^ is calculated using only the calorimeter and is not affected by 
the muon momentum resolution) and J^t > 20 GeV. An isolation condition, A7l{e, fi) > 0.5, 
where A7l{e, fi) is the separation between electron and muon in rj-cp space, is applied to 
remove background from cosmic ray muons accompanied by a bremsstrahlung photon. The 
/i/i candidate events are selected by requiring a leading muon with px > 25 GeV/c and a 
second muon with pt > 20 GeV/c. The projection of J^t onto the bisector of the muon 
tracks in the transverse plane is required to exceed 30 GeV in order to remove background 
from Z boson production. This is less sensitive to the muon momentum resolution than a 
dimuon invariant mass requirement. 

Additional cuts are applied similarly to all three channels. To remove background from 
Z TT, and from bb production in the fi/i channel, the transverse opening angle A0f 
between one charged lepton i and I^t is required to be less than 160°. This cut is applied 
to the second-leading electron in the ee channel, the muon in the e/i channel, and the 
leading muon in the /in channel. In addition, for the ee and e/i channels, Acp^ is required 
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to exceed 20°, and both A0£ requirements are removed if Ifx > 50 GeV. Also, to reduce 
background in all three channels from tt production, the hadronic Et in the event, E^'^'^ = 
— {E^ + E^ + I^t), is required to satisfy E^^^*^ < 40 GeV. After imposing these selection 
criteria, one ee candidate, two e/x candidates, and one fifi candidate remain. 

The detection efficiencies for W boson pair production with SM and anomalous cou- 
plings are determined using a fast Monte Carlo program (the Monte Carlo event generator 
of Ref. |T^ plus a parametric detector simulation). The detection efficiencies for SM W 
boson pair production are also calculated using the PYTHIA event generator followed 
by a detailed GEANT simulation of the D0 detector and are found to agree with those 
determined from the fast Monte Carlo. Trigger and particle identification efficiencies are 
determined from the data. The trigger efficiency for the ee and e/i data samples is (99l3)%. 
For the /i/i sample, the trigger efficiency is (68.7 ± 5.8)%. Table | shows the detection effi- 
ciencies for SM W boson pair production events and the number of expected events based on 
a cross section of 9.4 pb |T^. The systematic uncertainty in the detection efficiency comes 
from electron (2.2%) and muon (7.5%) identification, electron (2.0%) and muon (8.5%) trig- 
ger efficiencies, and the difference between the detection efficiencies estimated with the two 
Monte Carlo methods (5%). Uncertainty due to the choice of parton distribution function 
and evolution scale (5%) is included in the uncertainty on the number of expected events. 

Backgrounds due to Drell-Yan dileptons, W'-f, tt, and Z boson production are estimated 
using the PYTHIA, ISAJET p| and HERWIG p| Monte Carlo event generators, followed 
by the detailed GEANT simulation of the D0 detector. Backgrounds due to high-p^ Z ^ ee 
and /i/i events are studied using a Monte Carlo event generator based on the theoretical 
model of Ref. PO] and the parametric detector simulation. Backgrounds from multijet and 



W + jet events with a jet misidentified as an electron and with heavy quark production of 
isolated muons are estimated from the data. The probabilities for a jet to be misidentified 
as an electron and for a jet to be accompanied by a muon that satisfies the isolation criterion 
are measured from a large sample of events passing jet triggers. Events with large are 
rejected from this sample to remove ly+jets events. For electrons, the misidentification 
probability is found to be a slowly rising linear function of jet E^ (3.7 x 10~^ at 20 GeV, 
1.9 X 10-^ at 100 GeV in the central region; and 3.5 x 10"^ at 20 GeV, 1.8 x 10"^ at 100 
GeV in the forward region), while for muons it is found to be constant (1.5 x 10^^ for the 
e/i sample and 1.5 x 10~^ for the /i/i sample). The background estimates are summarized in 
Table ||. Systematic uncertainties include those listed above as well as the uncertainty on 
the production cross section of the background processes. 

The number of candidate events, four in the 1993-1995 data sample (five when the 1992- 
1993 data sample with one ee candidate and 0.6 ±0.1 background events is included), is 
consistent with an expected SM WW signal of 1.5 ±0.1 (1.9 ±0.1) events plus an estimated 
background of 2.5 ± 0.4 (3.1 ± 0.4) events. We have studied the stability of the results 
by relaxing some of the event selection requirements, e.g. E^^*^ and dielectron invariant 
mass criteria. The increases in the numbers of candidate and background events are found 
to be consistent with the expectations. An upper limit on the W boson pair production 
cross section is calculated from the number of the candidate events and the estimated back- 
ground events using the Poisson-distributed number of events convoluted with Gaussian 
uncertainties on the detection efficiencies, background, and luminosity. For SM W boson 
pair production, the upper hmit for the cross section is 37.1 pb at the 95% C.L. using the 
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1992-1993 and 1993-1995 data samples. The probability that the observed number of events 
correspond to a fluctuation of the background, with no signal, is 20.6%. 

By studying Et spectra of leptons from W boson pair candidates, limits can be obtained 
on the anomalous WW'y and WWZ couplings. Use of this kinematic information provides 
significantly tighter constraints on anomalous couplings than those from the measurement 
of the cross section (the method used in previous WW — > iFuu' analyses P,^]), since the 
predicted increase in the gauge boson pair production cross section with anomalous cou- 
plings is greater at higher gauge boson px- A binned maximum likelihood fit is performed 
to the measured spectra of Et of the two leptons in the event. Two-dimensional bins in 
Et of one lepton versus Et of the other lepton are used in order to take into account the 
correlation between the two leptons in the event. The binnings used in the fit are shown 
in Table |ITI|. The probability for the sum of the background estimate and Monte Carlo 
WW signal prediction to fluctuate to the observed number of events is calculated in each 
bin for a given set of anomalous couplings. The uncertainties on the background estimates, 
efficiencies, integrated luminosity, and theoretical prediction of the WW production cross 
section are convoluted with Gaussian distributions into the likelihood function. The likeli- 
hood functions are calculated for the 1992-1993 and 1993-1995 data samples separately and 
are combined taking into account correlated uncertainties, such as theoretical uncertainties. 

The WW production process involves the WW'j and WWZ couplings, unlike the W'j 
production process which depends only on the WWj couplings. Limits on anomalous cou- 
plings are obtained using two assumptions on the relationship between the WW'j and WWZ 
couplings. Figure shows bounds on anomalous couplings from this analysis and from the 
unitarity condition p!^j21[] using A = 1.5 TeV. In Fig. |l](a), the values for Ak and A are 



assumed to be equal for the WW'-) and WW Z couplings. Limits at the 95% C.L., when A 
or Aft is set to zero, are: 

-0.62 < Aft < 0.77 (A = 0); 
-0.53 < A < 0.56 (Aft = 0) 



In Fig. 0(b), the HISZ relations |jTT| are used. Limits at the 95% CL. using the HISZ 
relations are: 

-0.92 < Aft^ < 1.20 (A^ = 0); 
-0.53 < A^ < 0.56 (Aft^ = 0) 

The innermost curve is the 95% CL. contour when only one coupling is treated as a free 
parameter (e.g., limits on the axes) while the middle curve is the 95% CL. contour when 
both of the couplings are free parameters. All of the limits obtained in this analysis are 
comparable to those obtained from the analysis of WW/WZ — >■ eujj events 

In summary, a search for WW ii'uu' in pp collisions at y/s = 1.8 TeV is performed 
using the 1992-1993 and 1993-1995 data samples. In approximately 97 pb~^ of data, five 
candidate events are found with an estimated background of 3.1 ± 0.4 events. From the 
standard model, 1.9 ±0.1 events are expected. The number of observed events is consistent 
with the standard model prediction plus background estimate. The 95% CL. limits on the 
anomalous couphngs —0.62 < Aft < 0.77 (A = 0) and —0.53 < A < 0.56 (Aft = 0) are 
obtained from a binned maximum likelihood fit of the Et spectra of leptons, assuming equal 
WW'J and WWZ couplings. 
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FIG. 1. Contour limits on anomalous couplings for A = 1.5 TeV: (a) Ak = Ak^ = Akz, 
X = X-y = Xz', and (b) HISZ relations. The innermost and middle curves are 95% C.L. one- and 
two-degree of freedom exclusion contours from the fit of the Et spectra of leptons, respectively. 
The outermost curve is the constraint from the unitarity condition. Monte Carlo statistics limit 
the accuracy of the contours to ±0.01. 

We thank the staffs at Fermilab and collaborating institutions for their contributions 
to this work, and acknowledge support from the Department of Energy and National Sci- 
ence Foundation (U.S.A.), Commissariat a L'Energie Atomique (France), State Committee 
for Science and Technology and Ministry for Atomic Energy (Russia), CAPES and CNPq 
(Brazil), Departments of Atomic Energy and Science and Education (India), Colciencias 
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TABLES 



TABLE I. Detection efficiencies and SM signal event expectations for the 1993-1995 data 
sample. The uncertainties include both statistical and systematic contributions. 



Channel 


ee 


e/j, 




Detection Efficiency (%) 


6.03 ± 0.36 


4.76 ± 0.49 


1.19 ±0.18 


SM Expectation (events) 


0.52 ± 0.04 


0.86 ±0.10 


0.09 ± 0.01 


TABLE II. Summary of backgrounds and candidates for the 1993-1995 data sample. The units 


are number of events in the data sample. The uncertainties include both statistical and systematic 


contributions. 










ee 


e/x 




Background: 








Z — > ee or fi/j. 


0.27 ±0.06 




0.39 ± 0.09 


Z ^TT 


0.10 ±0.07 


0.21 ±0.08 


< 10-=^ 


Drell-Yan dileptons 


0.03 ± 0.04 




< 10-3 


Wj 


0.18 ±0.07 


0.35 ±0.14 




tt 


0.13 ±0.05 


0.18 ±0.06 


0.02 ± 0.01 


multijets/VF + jets 


0.20 ±0.14 


0.43 ± 0.28 


0.03 ± 0.01 


Total background 


0.91 ±0.19 


1.17 ±0.33 


0.44 ± 0.09 


Data 


1 


2 


1 



TABLE III. The binnings used in the maximum likelihood fit to set limits on the anoma- 
lous couplings and the numbers of candidate events (background estimate) for the 1992-1993 and 
1993-1995 data samples. 



ee channel (96.6 ± 4.5 pb ^) 


Tpel \ Tpe2 

H/rp \ tLlrp 


20 - 40 (GeV) 


40 - 500 (GeV) 


25 - 40 (GcV) 


2 (0.50 ±0.10) 




40 - 500 (Gc}V) 


(0.35 ± 0.07) 


(0.27 ± O.OG) 


e/x channel (96.2 ± 4.5 pb-^) 


\ EI^ 


15 - 40 (GeV) 


40 - 500 (GeV) 


25 - 50 (GcV) 


2 (0.95 ± 0.27) 


(0.16 ±0.05) 


50 - 500 (GeV) 


(0.16 ±0.05) 


(0.16 ±0.05) 


nil channel (77.4 ± 3.6 pb ^) 


Et^.' \ Ef 


20 - 40 (GeV) 


40 - 500 (GeV) 


25 - 40 (GeV) 


1 (0.08 ± 0.02) 




40 - 500 (GeV) 


(0.18 ± 0.04) 


(0.26 ± 0.05) 
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